gusucode.com > MATLAB与控制系统仿真实践 源代码 > MATLAB与控制系统仿真实践 源代码/第十四章/eg14_4.m
%绘制系统的nyquist图 num=[20 20 10]; den=conv([1 1 0],[1 10]); nyquist(num,den) %进行局部放大 num=[20 20 10]; den=conv([1 1 0],[1 10]); nyquist(num,den) %绘制系统Nyquist曲线 axis([-2 2 -5 5]) %得到局部的nyquist曲线 num=[20 20 10]; den=conv([1 1 0],[1 10]); w=0.1:0.1:100; %指定频率范围 [re,im]=nyquist(num,den,w); %返回Nyquist曲线参数 plot(re,im) %使用Nyquist曲线参数绘制曲线 axis([-2 2 -5 5]); %指定横纵坐标 grid; title('系统(20s^2+20s+10)/[(s^2+s)(s+10)]Nyquist图(\omega>0)','fontsize',12); xlabel('Real');ylabel('Imag')